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We exhibit, in the form of some identities in law, some connections between tilted stable subordi- 
nators, time-changed by independent Gamma processes and the occupation times of Bessel spiders, or 
their bridges. These identities in law are then explained thanks to excursion theory. 



1 Introduction 

1.1 Aim of this work 

The genesis of this paper is our interest in a class of subordinators (T t ,t > 0) which enjoy the 
following properties: (i) (T t ,t > 0) is of GGC (Generalized Gamma convolution) type. That is, it 
has no drift and its Levy measure is of the form 

(1.1) 8—E{e- xG ] 

x 

for G, a random variable such that E[(log(l/G)) + ] < oo, and 9 > 0; (ii) each marginal law of T t , 
for fixed t, can be described explicitly. Of course, this second demand is not mathematically very 
precise but we are seeking descriptions such as: 

T t ~X t Y u for fixed t 

where X t and Y t and have classical (e.g. beta, gamma) distributions. [Throughout, we shall write: 

U ~ V, to mean that U and V are identically distributed.] As we see above the law of the subordi- 

(ft 

nator (Tt), which can be written more precisely as (T t ' ) depends both on G and 9. However, we 
note that (T t {e ' G) ,t > 0) is equivalent (in law) to {T^ G) ,t > 0). 



Remark 1 . The general class of GGC random variables has been extensively studied by Thorin (1977) 
and Bondesson (1992). These are infinitely divisible random variables and hence one can naturally 
construct subordinators from them. We note that in general the Levy measure of such GGC subor- 
dinators may be defined as x _1 J Q e~ xy U(dy), where U is a sigma-finite measure. That is to say it 
does not necessarily correspond to, or arise from, a random variable. The case where U(pc) = oo, 
contains many important classes of random variables, including positive stable random variables of 
index < a < 1 and the GIG (Generalized Inverse Gaussian) class. See Bondesson (1992) for many 
properties and examples of GGC random variables. 
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1.2 First examples 

Here, we next give two examples of subordinators which satisfy the desired properties (i) and (ii) 
mentioned above. First let us take 6 = 1/2, then: 

Example 1 G:= G1/2 ~ 0(1/2 1/2), where /3( a 5) denotes a beta random variable, with parameters 
(a,b). 

Example 2 G := l/G 1/2 - l//3(i/ 2 ,i/2)- 

We denote these subordinators by E\ and £2, respectively (where here e stands for example). Then, 
as we shall see in section 3, there are the equations; 

(1.2) E[exp(-Aei(f))] = (VT+A-VX)* ' 



(VTTA + Vx) 

and Ei(t) ~ 7t/2//3(i/2,(i+t)/2)- 

(1.3) E[exp(-Ae 2 (t))] = 

(1 + VTTA) 

and ei{t) ~ 7t/2/ 3 ((i+t)/2,(i+t)/2)- 

Remark 2. Interestingly, the processes ex an d £ 2, have appeared previously. e\ is a special 
case of the random variables recently studied in Bertoin, Fujita, Roynette, and Yor (2006) and 
Fujita and Yor (2006). The fact that exit) ~ 7t/ 2 //5(i/2,(i+t)/2), as indicated in ()1.2p . is demon- 
strated in Roynette and Yor (2006) where connections to Feller (1966) are also noted. The fact that 
e 2 (i) ~ 7t/2/%i+i)/2,(i+t)/2)) as indicated in (|1.3[) . follows from a result of Cifarelli and Melilli (2000), 
concerning Dirichlet mean functionals, as discussed in James (2006). 

1.3 A general class of variables G 

The fact that /3(i/2,i/2) is the arcsine law, and that this law is that of the time spent in 1R + (or 1R_) 
by real valued Brownian motion, led us to look for generalizations of the above examples, where 
the variable G in (|1.1[) would be replaced by a variable distributed as the time spent in H + by a 
symmetrized Bessel process and even to look for further generalizations. As shown in Barlow, Pitman 
and Yor (1989), the time spent in IR + by a symmetrized Bessel process with dimension d — 2(1 — a), 
for < a < 1, is distributed as; 

T a 
T Q + T' a ' 

where T a and T' a are two independent standard stable (a) random variables, i.e. E[exp(— AT Q )] = 
exp(— X a ). More generally, we looked for a description of the subordinators satisfying (jl.ip with 



(1.4) G = 



where all the coefficients are assumed to be non-negative, and the (T5 , j — 1, 2, ... , N) are N iid 
standard stable (a) variables. 

It turns out that the subordinator T t associated with G given in (|1.4| may be obtained from 

N 

%:=^' V) =X>4 , '' ) («) 

i=l 
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where the (Tq '(u);u > 0) are constructed from N independent stable (a) subordinators, each of 
them being Esscher transformed with parameter (fj), that is, 

E[exp(-AT^)(w))] = exp(- U [(^ + A) Q - <]). 

In fact, (T t ) is constructed as i^€ lt , t > 0) where (7*, t > 0) is a standard Gamma process independent 
of (^u,u > 0). We call the process (^ 7t ,i > 0) a (generalized) positive Linnik process with index 
a, and parameters (fJb,y). The name we use comes from its link to the class of random variables 
expressible as L Q , t = ^^"Ta, defined for each fixed t > 0, where j t and T a are independent. 
Naturally these random variables are equivalent in distribution to a positive stable (a) subordinator 
time changed by a gamma process 74, for each fixed t. The random variables, "L a ,t, are known in 
the literature [see for instance Bondesson (1992, p. 38) and Devroye (1990, 1996)] and are sometimes 
called generalized positive Linnik random variables. These random variables have a host of interesting 
properties and distributional representations. Huillct (2000) discusses some results related to the 
subordinator. James (2006) provides an extensive recent study of L Q t and its Esscher transformed 
version which connects with the work of Barlow, Pitman and Yor (1989), and our present exposition. 
In this sense, L Q t and its Esscher transformed version arise as special cases of (^ 7t , t > 0). Moreover, 
C^7t>* — 0)' a l so contains an interesting subclass of the models discussed in the recent related work 
of Bertoin, Fujita, Roynette and Yor (2006). We shall show in Section 4 that we are able to describe 
precisely the one-dimensional marginals of the generalized positive Linnik process, (^t> lt ,t> 0), which 
therefore answers positively our demands (i) and (ii) made in (1.1). These one-dimensional marginals 
are closely related with the occupation times of the Bessel spiders living in webs constituted by N 
rays as described in Barlow, Pitman and Yor (1989). Once this was noticed, it became of interest to 
explain this coincidence. This is what we do in the last two sections of the paper. 

1.4 Organization of the paper 

We now describe the organization of the remainder of this paper: in Section 2, we recall notation and 
results about Bessel spiders and their occupation times of rays; in Section 3, we show that we can 
describe the one-dimensional marginals and the Levy measures of the generalized Linnik process, 
in Section 4, we compare the results obtained in Section 2 and 3; in the final Section 5, we give an 
explanation, based on excursion theory, for the close relationship between the distributions we have 
encountered in Sections 2 and 3. 



2 On the occupation times of rays by a Bessel spider 

In this section, we simply recall the definition of a Bessel spider, with index a E (0, 1), and whose 
state space is E, the union of N half lines, or rays, ii, . . . , In, originating from a common point which 
we denote as 0. In fact, these random processes are indexed by both a € (0, 1), and p = (p\, . . . ,Pn), 
a generic probability on {1,2,..., TV}. Now, we may introduce y(a,p) = (S t ,t > 0), a Bessel spider 
with index a, and probability p, which lives on E, as, informally, a process (S t ,t> 0) which behaves 
like a BES(a) process on each of the rays, Ji, . . . , In, and which, when coming back to 0, rides off in 
all directions at once, and chooses its ray with probability pi. Of course, this is only a heuristic 
description as the point is regular for itself (with respect to the Markov process (St)), but rigorous 
constructions have been provided by Barlow, Pitman and Yor (1989), with the help in particular of 
excursion theory. 

The case a = 1/2 is particularly interesting, as the (1/2, p) spider behaves like a Brownian 
motion in each of its rays.; We also note that for a = 1/2, and N = 2, p = (pi,p 2 ), <5 fi (l/2,p) may 
be identified with the skew Brownian motion of, say, parameter p\ (or p 2 \). 
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We also recall that \S t \ := dist(S' t ,0) is a BES(a) process (in particular, for a = 1/2, it is 
a reflecting Brownian motion), and (S t ,t > 0) admits a local time at 0, (L t ,t > 0) such that 
(|>!?t| 2a — Lt,t > 0) is a martingale; we refer the reader to Donati-Martin, Roynette, Vallois, and 
Yor (2005) for a discussion of different choices of this local time found in the literature. In the 
following theorem we describe the joint laws of : {A^ = J Q * dsI(S s G 7j); 1 < i < N-.L t }, for fixed 

i, as well as the sequence {a[^;i < N,£ t }, the same quantity, but now considered for the standard 
spider's bridge. 

Theorem 1 (1) There is the identity in law 

i n 1/a tW 1 
{A\, 1< i < N} = { P * , f« - ,i< N; rr—-}, 

L^j y j a L^j rj -*- a 

where E[(exp(-AT ( J ) )] = cxp(-A Q ) 

(2) The law of{a^;i < N,(^ a } is absolutely continuous w.r.t that of; {(A\,i < N)]L^ a } with 
density T(l - a)L\. This may be written as: V / > 0, Borel ,: [0, 1]^ xl + ^ M + , 

1/a ^ ] .. ^ ,r. 1 \ r(l + a) . 



E[/((oW, i < N);£\' a )] = E[/( P ' ° ; 1 < TV 



In the next section, we shall see, in a different context, similar quantities-in particular, the ratios: 
T^/E^tW appear. 

3 An identity about multivariate positive generalized Linnik 
processes 

In this section, we consider 

^ ) (u) = (T^\u);l<i<N),u>0, 

an K^— valued process, whose components are independent z/j-Esscher transforms of stable(a)- 
subordinators. Let, furthermore, (7t,i > 0) denote a standard gamma process independent of ^i"^. 
The next theorem offers a description of the joint law of (^i (7* ) , 7t ) , in terms of a sequence of iid 
stable-(a) variables, (T^,i < N), and an independent gamma (to) variable. 

Theorem 2 1) For any F : — > M + , and 5 : M + -> K+, Sore/, one has: 
(3.1) 



E[ W M(^)]-^E - V( , 7 , .. y .^, 7 , , (i/ . Ta) 
w/iere cr = an d v ■ T a = Y^j=i v j ■ ^ n particular for a = 1, i/ie formula simplifies 



as: 



(3.2) E[F((^")( 7t ))ff(7t)] = Y^y E 



F[(-^-)T a ) 9 ((-^-) )- 



v T a y - yv V- r a ' '(„. T a y 



2) IfJtC(dxi, . . . , d^Ar) denotes the Levy measure of the N -dimensional Levy process (*£a (7t)j t > 
0) then for any F : — > R + continuous, with compact support in (0, 00) : 

(jSf.F) = lim(-E[ J P(^)( 7t ))]) = a / -e-«E[F( u _2_)] 

tlU C Jq U V ■ J. a 
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Remark 3. Note that the last formula follows from 

ft 



E[F(X; 



dsLF(X s )}, 



where (X t ) is a Levy process, starting at 0, F a regular function, equal to on a neighborhood of 
0, and L is the infinitesimal generator of X. If X has no drift, and no Brownian component, then 



LF{0) = (Sf,F). 



Corollary 1 In the case a — 1, the variables v-ffa \ft) and ^^{jt) / (jt) 1 ^" we independent and 
satisfy: 



*(")/ 



(i) v r <\lt)=l a t 



(ii) 



(3.3) 



E 



H 



tXit) 



v7t 



l/c 



T(l + at) 

r(i + t) 



E 



1 



(y ■ T a ) 



for any positive measurable function H 



or equivalently: 



(3.4) 



^(7 t )=^(l), 



7t 



where on the RHS of [3.4]) , the law of the variable appearing there is the mixture, with respect 
to the law of (vjl^ a ) of the laws o/ (1), i.e. the fi-Esscher transform o/^i ^(l). 7t is 
assumed to be independent of the process: (^i (l),/x > 0). 

Corollary 2 Let /Ltj > 0, > 0; 1 < i < N. The generalized positive Linnik process 

N 



is a GGC subordinator with Levy measure: 



" 1 E 



exp — x 



il 1 ' T) 



dx. 



Corollary 2 follows immediately from the expression of Jz? (or rather (Jzf , F)) in point 2) of Theo- 
rem 2. Likewise, Corollary 1 is easily deduced from formula (|3.2[) in theorem 2, by easy changes of 
variables. We now give a proof of theorem 2: 

Proof of Theorem 2: Of course, (|3.2p follows immediately from (|3.1[) . so it suffices to prove (|3.1[) . 
We first use the absolute continuity relationship between the laws of (^i (s) ; s < it) and (s);s < 
u), which yields that the LHS of (|3.ip equals: 

N 

E[F(Cri°)( 7t ), 7 t) exp{- E(^ 0) (7t) - 7t< }] = E[F( 7f 1/Q T Q , 7t ) exp(- 7t 1/c > • T Q ) + a 7t )]. 

»=i 

(by scaling and independence), 

= F7T / fa f " 1 e l(<T_1) E:[f(x 1/o T a ,j;)exp(-3: 1 / a ( i /.T Q ))] 

r(t) Jo 
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(then, changing variables x = y a ), 

/>oo 

= W) J d ^ Q " 1 (^) t_le " Q( ^ 1)E ^(2/ T - v a ) e M-y(" ■ T ))] 



(then changing variables: y — zj(v ■ T Q )), 
aT(pd) 



T(t) 

which yields the desired result. 



Remark 4 FKe afao no£e that the preceding argument yields, but we do not give the details, the 
following 1-parameter extension of formula \3.1\l 
(3.5) 



m 11 



r(l + at) 



E 



7at 



m ^( ( ^% )TQ ) 5(m( 



7at N a x eX P((li?3t) O ( '- m )) 



3.1 Examples 

We now illustrate Theorem 2 and Corollary 2 in the particular case N = 2 and a = 1/2, for which 
the law of T\/ 2 is explicit (recall E[e~ ATl / 2 ] = e _v/A ). We shall also assume yfv\ + = 1. 

Example 1. Set v x = 1, v 2 = 0, Mi = A*2 = 1; E[exp(-AT 4 )] = (vT+A - \/A) . 
Proposition 3.1 77ie Lew?/ measure of the process (T t ) is 

1 du 

2~ E i eX P(- W P(l/2,l/2))J 



fn) TTie distribution of the variable T t is that of 



7t/2 



Example 2. Set i/ x = = 1/4, ^ = 1/4, /i 2 = 0; E[exp(-AT t )] = (2/(1 + VT+A)) 
Proposition 3.2 (i) The Levy measure of the process (T t ) is 

1 du„ 



2 u 



•E[exp(-u//3 (1/2il/2) )] 



(ii) The distribution of the variable T t is that of 



7t/2/ 3 (i±*,i±i)- 



Example 3 (This is the general case for a = 1/2, N = 2): Set yjv\ + ^Jv^ = 1, fx 2 > 0; 
Proposition 3.3 (i) The Levy measure of the process (Tt) is 

1 du_ r . v\ Tx/2 + v 2 T x / 2 
-— Eexp-u 1 -J-))] 

2 u mT 1/2 + n 2 T 1/2 
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(ii) The law of the variable T t satisfies: V/ > 

Mi7(i+t)/ 2 + M27(i+i)/2' 



(3.6) 



E[/(T t )] = C t E 



/ 7t/2 



^7(l +t )/2 + »*7(l+t)/2 J (V!f a+t) +^27<1±£))* /2 



for some universal constant C t , and"/, 7, 7' are independent gamma processes. 

Rather than proving Propositions 3.1, 3.2, 3.3 in that order, we shall first prove Proposition 3.3, and 
then see that the results of Proposition 3.1 and 3.2 follow. In fact, we only need to prove the point 
(ii) in each of the Propositions. 
Proofs of points (ii): 

(a) From Theorem 2, we know: 



nm)} = 



r(i + t) 



/ It/: 



/2 



,. T (2) 
^2 1 1/2 



(•*t$, 



,. T (2) x*/2 



Since T 



(0 A 
1/2 _ 



c /li/2 f° r * — 1) 2, we get, with obvious notation: 



E[/(T t )] = C t E 



/ 7t/: 



7i7i 



MlTl/2 +M27l/2\ / 
' ^i7i/2 + ^271/2 y \(^l7l + ^27i) ; 



t/2' 



which simplifies to (|3.6|) . 

(b) In the case of Proposition 3.2, we have: v\ = and then 7(i +f )/ 2 /(7(i +t )/ 2 + 7(i+t)/2) is a 
BETA((l+t)/2,(l+t)/2) random variable, independent from (j[ 1+t y 2 +7(i+t)/2)- 

(c) In the case of Proposition 3.1, the formula (|3.6p in point (ii), Proposition 3.3., simplifies as: 



nm)] = c t 



f it/ 



E 



/ It/: 



%l+t)/2 + 7(l+t)/2 
^(l+t)/2 

7i/ 2 +7(l+t)/2\ 



(7'1-M) ) 



t/2 



= E 



/ 7t/s 



7i/ 2 
1 



3,1 i+t 

V 2 ' 2 



It might also be of interest to discuss the case N — 2 and general < a < 1 in the same vein, 
because the law of (T a /T' Q ) is simple, as shown by Lamperti (1958). Let us denote by R Q this ratio. 
Then we have the following: 

Proposition 3.4 (N = 2, < a < 1, vf + v% = 1 ). 

(i) The Levy measure of the process is: 



du w 
a — E 

u 



(ii) The law of T t satisfies V/ > 0, 



(3.7) 



E[/(T t )] - CtE 



I It/: 



viR a + v 2 
exp [ — u — — 

Ml-Ra + M2 
fJ,lR a + fJ>2 



1 



1 



where, on the RHS, 7 t / 2 is independent from the pair (T a , r t a ) and R a = T a / l! c 
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3.2 Looking for a simplification of Proposition 3.4 

Using a conditional expectation argument we can replace (T^)~ Qt appearing in (|3.7|) with 



h at (R a ) :=E[(T'J Qt |R Q ]. 

Naturally this expression has utility only if it has a tractable form. Here we obtain such an expression 
by using Kanter's (1975) explicit representation of the density of a positive stable random variable, 
in conjunction with Lamperti's expression for the density of R a . This form of the stable density is 
apparently not well-known. First let f a denote the density of a unilateral stable- (a) density. Then 
generically the conditional density of T^|R Q = r is given by 

, /is sf a (rs)f a (s) 

Now setting 

K a (u) 



sin(7rait)\ / sin((l — a)iru) 



sin(7ru) J \ sin(7rcra) 
it follows from Kanter (1975) [see also Devroye (1996)] that 



/q(s) = ^_ a -i/d-°0 / e' s ^ K ^K a (u)du. 
1 — a 



That is T' a = (K a (U)/e) {1 - a)/a , where U is a Uniform[0, 1] random variable independent of e, 
which is exponential(l). Additionally from Lamperti (1958)[see, also Chaumont and Yor (2003, p. 
116)] we obtain, 

sin(7ra) 

7r r 2a + 2r a cos(7ra) + 1 
Combining these points we arrive at the following result 



fn (r) = —z ; : for r > 0. 



Proposition 3.5 (i) The conditional density of l* a \R a — r is expressible as 

«(2-o0 2 



sin(7ra) \ 1 — a 
where C a (ui,u 2 \r) — [r~T=z K a {ux) + K a (u 2 )] 
(ii) It follows that h at (R a ) :=E[(I y Q ) Q *|_R Q ] equals, 

(3.8) h at (R a ) = — T^T : f ) Da.at(Ra) 

sm(7ra) \ 1 — a J 

where D a ^ at (R a ) = J* J* [C a (ui, u 2 \R a )]~ {t(1 ~ a)+2) K a (ui)K a {u 2 )duidu 2 . 

This leads to an alternative characterization of (Tt) in proposition 3.4 as follows; 

Proposition 3.6 The distribution of the variable T t in Proposition 3.4 with N = 2, < a < 1, 
z/f + u% = 1, satisfies V/ > 0, 



E[/(T t )] = C t E 
where h at (R a ) is given by (3.8\) . 



. /j,iR a + iJ, 2 \ h at (R a ) 

f 7t/2" 



ViR a + V 2 ) {viR a + V 2 ) 
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4 A comparison of Theorems 1 and 2 

This comparison is only partial, as here we shall use the hypothesis: a = Yli=i v t — 1) 80 tnat we 
can associate with the Esscher sequence {vi)i<N the probabilities {pi)i<N defined by pi — vf. Then 
a simple cross- inspection of Theorems 1 and 2 shows that, taking t = 1, 

(4.1) {(^)(e);i < iV),e} = {( 7 a4 <) ;< < (7a) a 4}, 

where, on the LHS, e is an independent exp(l) variable, and on the RHS, 7 Q is independent of the 
Spider's Bridge. We shall now transform the identity in law (|4.1[) . or rather its RHS, by using the 
two following facts: 

(i) if g = sup{s < 1 : S 3 = 0}, then g = /3 (q ,i_ q) , 

(ii) 7a = P(a,i-a) e where /9( Qi i_ Q ) denotes a beta variable with parameters (a, 1 — a), which is 
independent from e. 

Thus the identity in law (|4.ip may be rephrased as: 

(4.2) {fa«£* } (e);i < A0,e} = {(A«5* < iV),^}, 
where f/t = sup{s < t : S s = 0}, which is equivalent to: 



f (i) L e = e 

(43) I (it) {(4^;i<7V)|L e -/} = {(^i ,yi) (e);i<7V)} 

In the next section, we give an explanation of 



5 An explanation of (4.3) 

We shall now show that (4.3) follows from the conjunction of the next lemma and again the scaling 
property of a Bessel spider. 



5.1 Lemma 

Lemma 5.1 Let (X t ,t > 0) denote a nice Markov proces taking values in a general space E; let 
belong to E, such that is regular for itself (with respect to X), and denote by (L t ,t > 0) a choice 
of its local time at 0; let (t£,£ > 0) denote the inverse local time, i.e. 

T t = inf{t : L t > I}. 

It is a subordinator with associated Bernstein function, 4>(Q): 

E{e- 0T <] = e-^W. 

Then, if eg denotes an exponential (9) random variable independent from (X t ,t > 0), and (A tl t > 0) 
is a continuous additive functional, one has 

(i) L eg is exponentially distributed with parameter ij){ff); 

(ii) The conditional formula holds; 



\L ee =£} = 



E[e- 9T '-} 
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This lemma is quite classical (it may be obtained from excursion theory), let us admit it for a 
moment and derive (4.3) from it. For this purpose, we take: A t — S^Li XjA^\ Xj — 0, and we write 
(ii) in the lemma as follows: 

E[exp(- ]T XjAg] )\L eg = I] = E[exp(- £(Aj + 0)A%)] exp(W Q ) 

3 3 

(since: n = £\ 

= ex P (-iJ2 Pj ((\ j + er-e<*) 

(since the Ito excursion measure n associated with the local time (L t ) satisfies n = J2jPj n j> where 
rij is n restricted to excursions on Ij) 



E 



exp -^A^ffl 



and it now suffices to take 8 = 1. 



5.2 A global view of the lemma 

In fact, the lemma is only a part of the more general statement. 

Proposition 5.1 (We use the same notation as in the Lemma). The following holds: 

(i) The pre-g ee process {X u ,u < g ee } and the post-g eg process {X geg+V ,v < d eg — g ee } where 
d t — inf{s > t : X s — 0}, are independent; 

(ii) The law of (X u , u < g eg ) may be described as follows: 

(a) L ee = Lg se is exp(i/>(6*)) distributed; 

(b) {(X u ,u < g eg )\L eg = 1} is distributed as: (X u ,u < Tg) under the probability 

exp(-6»r £ + lijj(6)) ■ P 

(Hi) The law of the process {X 9e +v , v < d eg — g eg } is 

1 



-(1 -exp(-0V(e)))n(de) 



where V(e) denotes the lifetime of the generic excursion e. 

For details of the proof and further references, the reader may consult Salmincn, Vallois and 
Yor (2006). 
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